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The multivariate central limit theorems (CLT) for the volumes of excursion sets of stationary 
quasi-associated random fields on R d are proved. Special attention is paid to Gaussian and 
shot noise fields. Formulae for the covariance matrix of the limiting distribution are provided. 
A statistical version of the CLT is considered as well. Some numerical results are also discussed. 
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An important research domain of modern probability theory is the investigation of ge- 
ometric characteristics of random surfaces (see, e.g., [1-3]). The origin of interest often 
roots not only in pure mathematical challenges but also in various applications, including 
those in industry. Wc mention one motivating example for our study. 

The contemporary method of papcrmaking goes back to the Han Dynasty period. 
Nowadays, the method is essentially the same, but machines in modern pulp and paper 
mills operate much faster. The surface structure of the paper during the forming process 
S^ ■ determines the quality of the production. 

H , To model the paper surface, stationary random fields, say, shot noise (cf. [4]) or Gaus- 

sian, can be a reasonable first choice. Comparing by eye real paper image data and 
simulated realizations of such fields, one easily concludes that the similarities are strik- 
ing. But it is hard to quantify how different these two images really are. To test whether 
the available image data originate from a realization of a specified stationary random 
field, the excursion sets can be considered. 

We prove the central limit theorem (CLT) for volumes of excursion sets of a stationary 
field X = {X(t),t£ T}, T C M d , to characterize the surface generated by X. It is reason- 
able to assume that the field X could possess a dependence structure more general than 
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positive or negative association used in a number of stochastic models; see, for exam- 
ple, [7]. Our main results yield uni- and multivariate CLT for quasi-associated random 
fields. The CLT is generalized in [12], page 80, having been obtained by other methods 
for volumes of excursion sets of stationary and isotropic Gaussian random fields. We also 
discuss the consistent estimators for the asymptotic covariance matrix that arises in the 
limiting distribution. 

Note that we do not tackle here the interesting problems concerning the study of 
moving levels for excursion sets, the estimate of the convergence rate to the limit law 
and the analysis of the functionals in Gaussian random fields based on the Dobrushin 
Major techniques. In this regard, we refer, to [12, 15-17]. 

As to the problem of characterizing the paper quality taking into account the "hills" 
and "valleys" of its surface discernible with the help of microscope, it is by no means 
simple. In fact, we have to specify the admissible (average) number of such hills along 
with their size. Moreover, the thickness of the paper should be controlled as well (no 
holes or high peaks) . Thus the study of the excursion sets for random fields is the first 
natural step to investigate such random surfaces. The application to paper surface image 
data will appear in a separate paper. 

The present paper is organized as follows. Section 2 provides preliminaries on depen- 
dence concepts related to association and excursion sets for random fields. The CLT 
for the volumes of excursions of quasi-associated stationary random fields over one or 
finitely many levels are formulated and proved in Section 3. The special cases of station- 
ary shot noise and Gaussian random fields are treated in more detail. Section 4 contains 
a statistical version of the limit theorems mentioned above where the (unknown) limit- 
ing covariance matrix is consistently estimated. Numerical results illustrating the limit 
theorem of Section 3 are given in Section 5. Finally, we conclude with the discussion of 
some open problems. 



2. Preliminaries 

In this section, we recall some dependence concepts for systems of random variables. 
Various examples can be found in [7]. After that, we introduce the excursion sets that 
are the main objects of this study. Then we consider the sequences of regular growing 
sets forming observation windows. 



2.1. Dependence concepts for random fields 

Consider a family, X = {X(t),t e T}, of real- valued random variables, X(t), defined 
on a probability space, (fi, T, P). A set, T, will be a subset of R d or Z d . For I C T 
let Xj = {X(t),tE I}. Introduce the class M(n) consisting of real-valued, bounded, 
coordinate- wise non-decreasing Borel functions on W l ,n <E N. The cardinality of a finite 
set, U, will be denoted by card U. 
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Definition 1. A real-valued random field X = {X{t),t £ T} is called positively asso- 
ciated (we write X £ PA) if, for every disjoint finite set I,JcT and any functions 
f £ M(cardl) and g £ M(cardJ), one has 

cov(/(X / ), 5 (X, / ))>0. (1) 

Here, we use any permutation of (coordinates of) the column vector (X(ti),. . . ,X(t n )) T 
for Xi, I = {ti, . . . , t n } C T (and the analogous notation is employed for Xj ) ; T stands for 
transposition. Definition 1, given for any (not necessarily disjoint) subsets / and J C T, 
introduces the family of associated random variables (X £ A). The change of the sign of 
inequality (1) leads to the definition of negative association (one writes X £ NA). Clearly, 
X £ A implies X £ PA. Any collection of independent random variables is automatically 
PA and NA. Due to Pitt [19], a Gaussian family X = {X(t),t £ T} of random variables is 
associated if and only ifcov(X(s),X(t)) > for all s,t £ T. For such families, the concepts 
of A and PA coincide. A theorem by Joag-Dev and Proschan [13] states that a Gaussian 
family X = {X(t),t £T} £ NA if and only if cov(X(s),X(t)) < for s,t£T, s^t. 

Let BL(n) denote the class of bounded Lipschitz functions / : R™ — > R (n £ N) and 

Lip(/) = sup l/ i | X) ~y <oo, ||z||i = Y> fc |, x = ( Xll ...,x n ) T £« n . 

x¥:y \\x-y\\i £-J 

Since all norms are equivalent in R™, we sometimes use the Euclidean norm ||x||2 = 
(^22=i x "k) 1 ana - ^ ne supremum norm ||x||oo = naaxfc = i ] ... ]TC \xu\ of x £ R™ for the sake of 
convenience. 

Definition 2. A random field X = {X(t),t £ T} consisting of random variables X(t) 
with EX(t) 2 < oo is called quasi-associated (X £ QA) if 

|c v(/(X / ), 5 (X J ))|<Lip(/)Lip( g )^^|cov(X(.s),X(t))| (2) 

sei teJ 

for all disjoint finite sets I,JGT and any Lipschitz functions 

/:R card/ ^R and ,g:R cardJ ^R. 

If X £ PA or X £ NA and EX(t) 2 < oo for all t £ T, then (2) holds as was proved 
in [9]. Every Gaussian random field X (with covariance function taking both positive 
and negative values) is quasi-associated; see [20] and references therein. 

Definition 3. A real-valued random field X = {X(t),t £ Z d } is called (BL, 6) -dependent 
(X £ (BL,8)) if there exists a non-increasing sequence 6 = (# r ) r gN,0r —* as r — > oo, 
such that, for any finite disjoint sets I, J C 1 d with dist(7, J) = r and any functions 
f £ BL(card I) , g £ BL(cardJ), one has 

| cm(f(Xi),g(Xj))\ < Lip(/) Lip(g)(card I A card J)0 r , (3) 

where dist(7, J) = min{||s — t\\oo: s £ I ,t £ J}. 
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If X = {X(t),t£ Z d } e QA, then X e (BL, 9) whenever the Cox-Grimmett coefficient 



u r 



sup J2 \cov(X.,X t )\ (4) 



s£L 't: ||s-t||oo>r 



tends to zero as r — > oo. In this case, one can take 6 r = u r in (3). 

For a random field X = {X(t),t £ K d } we use (see [5]) the following extension of (3). 
Let T(A) = {(ji/A, . . . jd/A): (ji. • • • , Jd) e 2 d }, where A > 0. 

Definition 4- A real-valued random field X — {X{t),t £R d } zs called (BL, 8) -dependent 
if there exists a non-increasing function 9 = (9 r ) r> Q,9 r — > as r — > oo, suc/i i/iai, /or 
a// A large enough and any finite disjoint sets I , J C T(A) tuii/i dist(/, J) = r, and any 
functions f <G BL(card 7), g <G BL(card J), one has 

|cov(/(X 7 ), 3 pO))| < Lip(/)Lip( 5 )(card/Acard J)A d 9 r . (5) 

In many cases, one can use the integral analog of (4) for 9 r . Thus for a (wide-sense) sta- 
tionary random field X — {X(t),t E M. d } G QA, having covariance function R(t), t £ R d , 
absolutely directly integrable in the Ricmann sense (i.e., when d=l; see, e.g., Feller [11], 
page 362. For d > 1, the definition is quite similar. One takes the partition of R d gener- 
ated by partitions of each coordinate axis and forms the corresponding upper and lower 
Ricmann sums.), relation (5) holds with 

r = 2 [ \R(t)\dt, r>0; (6) 

•/|N|oo>r 

see [5]. We shall also write 9(X) = 9 r (X) to emphasize that 9 in (3) or (5) refers to the 
field X. 



2.2. Excursion sets 

Now we recall the definition of an excursion set and illustrate it by Figure 1. 

For a real- valued random field X = {X(t),t £M. d }, we assume the measurability of X(-) 
as a function on M. d x ft endowed with the cr-algebra B(M. d ) ® T . 

Definition 5. Let X be a measurable real-valued function on R d and T C M. d be 
a (Lebesgue) measurable subset. Then, for each u £ M, 

A u (X,T) = {t£T: X{t)>u] 

is called the excursion set of X in T over the level u. 

Let Vd{B) be the volume (i.e., the Lebesgue measure) of a measurable set B C K d 
and I{C} denote the indicator of a set C. 
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Figure 1. Realization of a stationary centered Gaussian random field X with covariance func- 
tion co\/(X(0),X(t)) = exp(— 1| i || 2 ) (left figure), bright colours indicate high values of X. The 
excursion set A U (X,T) for u = is shown in black (right figure). 

Since X is measurable, the volume of the excursion set 

v d (A tt (X,T))= f l{X{t)>u}dt 



is a random variable for each «eK and any measurable set T C M. d . 

2.3. Growing sets 

Denote the boundary of a set B C M. d by dB. The Minkowski sum of two sets, A, B C K d , 
is given by A® B = {x + y: x G A,y G B}. The following concept of "regular growth" for 
a family of subsets in M. d will be used in the sequel. 

Definition 6. A sequence, (W n ) n( =fi, of bounded measurable sets, W n C M. d , tends to 
infinity in the Van Hove sense ( VH-growing) if, for any e > 0, one has 



mr\ j v d (dW n ®B e (0)) 

v d {W n )^oo and — ^ , ^0 

l>d{W n ) 



(7) 



as n — > co, where B e (0) ={i£ K d : ||x||2 < £} is the closed ball in M. d with center at the 
origin G M. d and radius e. 

UW n = (a(n),b(n)] = (ai(n),&i(n)] x • • -x (ad(n),bd(n)} is a parallelepiped, then W n —> 
oo in the Van Hove sense if and only if bk{n) — o,k{n) — > co as n — ¥ co for k = 1, . . . ,d. 

Definition 7. A sequence of finite sets U n C Z d tends to infinity in a regular way if 

cardSU, 



card U n 







as n — > co; 



(8) 



c/. C7;. ffere SU n = {j G Z d \ U n : dist(j, C/ n ) = 1} and dist(j, [/„) = min feeC/ „ ||j - fc|| 
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3. Central limit theorem 

Now we state and prove a CLT for the volume of excursion sets of random fields. Ivanov 
and Leonenko [12] studied stationary and isotropic Gaussian random fields. In our ap- 
proach, the isotropy of Gaussian fields is not required. Moreover, we consider a more gen- 
eral class of random fields possessing the quasi-association property. To avoid long formu- 
lations, we introduce the following two conditions for a random field X = {X(t),t £ K d }. 

(A) X is quasi-associated and strictly stationary such that X(0) has a bounded den- 
sity. Assume that the covariance function of X is continuous and there exists some a > 3d 
such that 

|cov(X(0),X(t))|=O(||£|| 2 - Q ) as||t||2->oo. (9) 

(B) X is Gaussian and stationary. Suppose that its continuous covariance function 
satisfies (9) for some a> d. 

Notice that continuity of the covariance function of X implies the existence of a mea- 
surable modification of this field. We consider only such versions of X. We exclude the 
trivial case when X(t) = const a.s. for all t £M. d . 

3.1. Quasi-associated random fields 

To prove the CLT for the volume of excursion sets of a random field satisfying condi- 
tion (A), we need the following auxiliary result. 

Lemma 1 ([7], Lemma 7.3.4). Let {U, V} £ QA, where random variables U and V 
are square-integrable and have densities bounded by a > 0. Then 

| cov(I{[/ > u},l{V > v})\ < 3 • 2 2/ V /3 | cov([/,l/)| 1/3 

for arbitrary u, v £ R. 

Theorem 1. Let X = {X(t),t £ R d } be a random field satisfying condition (A). Then, 
for any sequence of VH-growing sets W n C K and each MgK, one has 

vjiMXM) - Vd {w n nx(o)>u) i„ j n _^ 

\Jv d {W n ) 

Here — ¥ denotes convergence in distribution, Y u being a Gaussian random variable with 
mean zero and variance 

<7 2 (u)= f cov(I{X(0)>u},I{X(t)>u})dt£R + . (11) 

Proof. Fix any u £ K and transform a random field X into a field Z = {Z(j),j £ Z d }, 
setting 

Z(j)= [ I{X(t)>u}dt-P(X(0)>u), (12) 
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where the unit cubes 

Qi = {x = (ari, . . . , x d ) T £ R d : < x k < 1, k = 1, . . . , d} © {j}, j = (j u . . .,j d ) T £ Z d . 

The Fubini theorem implies EZ(j) = for any j £ Z d . It is easily seen that the field Z is 
strictly stationary and square-intcgrablc. Introduce 

J- = {] £ Z d : Q j C W n }, J+ = {j £ Z d : Q j C\W n ^0} (13) 

and 

W~= |J Qj, W+= |J Qj. 

jeJn je./+ 

Due to (7), we conclude (see [7], Lemma 3.1.2) that 

v d(W~)-> oo and Vd(W~)/va(Wn)-*l as n -> oo. (14) 

Write 

v d {A u {X, W n )) - u d (W n )P(X(0) > u) 



VMWn) 

v d {A u (X,W~)) - v d {W-)P{X{Q) > u) 

\ftdWn) 

, u d (A u (X,W n )) - v d (A u (X,W-)) - (v d (W n ) - u d (W-))P(X(0) > u) 



(15) 



We prove that the second term on the right-hand side in (15) tends to zero in probability. 
By Chebyshev's inequality, it suffices to show that 

var(v d (A u (X, W n )) - v d (A u (X, W~)))/v d (W n ) -► 0, n ->■ oo. 

Set 

Y n (j) = ! l{X(t) >u}dt- UdiQj n W n )P(X{0) > u) for j £ Z d , n £ N. 

Note that Y n (j) = Z(j) for j £ J~ and n £ N (clearly F n (j) and Z(j) depend on it as 
well). Applying the Fubini theorem and Lemma 1, we get 

var(u d (A u (X, W n )) - v d {A u {X, W~))) 

= vaiY J2 Y ^) 

je,J+\J- 

< V / |cov(I{X(s)>w},I{X(t)>u})|dsdi (16) 

_l. — JQ,xQ„, 



j,mGJ+\J, 
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<MW+\W~)J2 f GilcovfX^.Jr^l^dsdi 

<v d (W+\W-)(c 2 +C 3 JT J2 I \\s-t\\- a/ ^sAt\ 

\ r=r oj£ 1 d : \\j \\ x =r J Qo XQ i ) 

< u d (W+ \ W~) \C 2 + CiJTr^r-*/ 3 ) = C 5 u d (W+ \ W~) 

for some r$ > and all n £ N. The factors Ci do not depend on n. We used the inequality 
| cov(X(s),X(t))| < r 2 for all s,t £ R d , which is satisfied as varX(t) = t 2 for any t £ R d . 
We also took into account that 



d-l 



card{j€Z d : \\j\\oo = r} <C 6 r 

for each rgff and employed the inequality (9) with a > 3d. 

By (14), (16) and in view of the relation v d (W~) < v d (W n ) < v d {W+), we get 



va<r((v d (A u {X, W n )) - v d {A u (X, W-)))/y/v d (W n )) -»• 0, n -»• oo. 

Now we show that J~ , introduced in (13), tends to infinity in a regular way as n — > oo. 
Indeed, J~ C J n C J+, where J„ := W n r\Z d , n £ Z d . Due to [7], Lemma 3.1.5, J„ tends to 
infinity in a regular way. Thus, it suffices to mention that S J~ C SJ n U ( J n \ J~ ) and apply 
the relations card <5 J n / card J„ — > and card J+/ card J~ — >• 1 as n — > oo. Lemma 3.1.6 of [7] 
implies that W,7 = Uj£j _ Oj tends to infinity in the Van Hove sense as n —$■ oo. 

So, while establishing (10), we can assume w.l.g. that W n = W~ , that is, W n is a finite 
union of cubes Qj (n £ N) and the sequence (W n )neN is VH-growing. 

Observe that 

^ W n )) - ^(W»)P(X(0) > u) = E jeWn n^ Z U) , = ^ 
V^d(W„) v/i/# n ) 

As X = {X(t), t € M d } e QA, it follows from (6) and (9) that X £ (BL, 9) with 

6» I .(X) = 0(r- a+d ) asr->oo(r>0). 

For 7 > (and u fixed) introduce the Lipschitz functions /i 7 : R — > R by the formula 

{0, if x < u — 7, 

(x — u + -f)/j, if u — 7 < x < u, (17) 

1, otherwise. 

Superposition of two Lipschitz functions is also a Lipschitz one. Thus, for n £ N and 
7 > 0, the random field Z na = {Z nn (j),j = (ji, . . . , j d ) T £ Z d } £ (BL, 9), where 

Wj) = ^ E ^(^(ii + -,...,^ + -))-EM^(o)) ( 18 ) 

fe 1 ,...,fe d =i ^ ^ n n 'J 
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and the terms of a sequence 9(Z n _ 7 ) admit the estimate 

8 r (Z nn )<C 7 ^ 2 r~ a+d , reN, (19) 

with CV depending neither on 7 nor on n. 

It is not difficult to verify that the finite-dimensional distributions of the fields Z nn 
weakly converge to the corresponding ones of the field Z 7 as n — > 00, where 



Z 1 {j)= h^X(t))dt-Eh^(X(0)), jeZ d . (20) 

Consequently (see [7], Lemma 1.5.16), we can claim that Z 1 <G (BL,8) and 9 r (Z 1 ) is 
bounded by the right-hand side of inequality (19). Theorem 3.1.12 of [7], guarantees 
that, for each 7 > 0, 

■Sn(7):= , 7==. ; ^K u , 7 -AA(0,a 2 (M,7)), n->oo, (21) 

\Jv d (W n ) 

where 

cr 2 (u, 7 )= V cov(Z 7 (0),Z r (j'))= / cov(/i 7 (X(0)),/i 7 (X(t)))dteR+. 

Therefore, to prove (10), two steps remain. First of all, we estimate the difference of the 
characteristic functions of the random variables S n (j) and S n and show that it tends to 
zero as7^ 0+. After that, we verify that 

cr 2 (u,7)^cr 2 (u) as 7^0 + . (22) 

Set h(x) = I{x > u] and H 1 {x) = h 7 (x) — h(x), where x £ K. (and u G R is fixed). Then, 
for each A € R, one has 

| Ee iAs„ (7 ) _ Ec iA5„! < | A |e| 5b(7) _ Sn \ < |A| (^^y) ' , (23) 

where i 2 = — 1 and 

K(7) = e( 2 / (^(^W)-E^ 7 (A(t)))dt] . 

It is easily seen that 



V n (j)<v d (W n ) f \cm(H^{X{0)),H 7 {X{t)))\ 



dt. (24) 



Furthermore, we have 

|cov( J ff 7 (X(0)),i/ 7 (X(t)))|<(E(i/ 7 (X(0))) 2 E(i/ 7 (X(t))) 2 ) 1 /2< a7 , 
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where a is a constant that bounds the density of -^(0). If | co\/(X(0),X(t))\ 1 ^ 3 < 7, then 
reasoning similar to that proving Lemma 1 leads to the inequality 

|cov(ff 7 (X(0)),i/ 7 (X(i)))|<C(a)|cov(X(0),X(i))r /3 (25) 

with some C(a) > 0. Write a = 3(d + (m), h> 0, and take R = C7~ 1 /( d+ ' i ), where c> 0. 
Then, in view of (9) and due to the appropriate choice of c, one can conclude that for 
all 7 > small enough, 



F( 7 ):= / |cov(iJ 7 (X(0)),F 7 (X(t)))|d< 

V (26) 

< ai uj d R d + C s [ ||i||r /3 di < C,Y /{d+ ^\ 



\\t\\ 2 >R 

where uJd = n d ^ 2 /(T(d/2 + 1)) is the volume of the unit ball in M d with the Euclidean 
norm. Consequently, inequalities (23), (24) and (26) imply that the laws of 5 r „(7) and S n 
are close for all n large enough if 7 > is small enough. 

Next, we proceed to (22). By the arguments leading to (26) and invoking Lemma 1, we 
deduce that a 2 (u) < 00. Similar to (25), one shows that if | cov(X(s),X(t))\ 1 / 3 < 7, then 

|cov(/ l (X( S )), J ff 7 (X(t)))|<^( a )|cov(X( S ),X(i))| 1 /3 (27) 

with D(a) > depending on a only. The absolute value of a 2 (u,y) — <J 2 {u) does not 
exceed the following expression: 

F(i)+ f \cov(h(X(0)),H J (X(t)))\dt+ f \cov(H^(X(0)),h(X(t)))\dt. 

Taking into account the above upper bound and relations (26) and (27), we complete the 
proof of (22). The asymptotic (finite) variances a (11,7) are non- negative, whence one 
concludes that cr 2 {u) > 0. 

In view of (21)-(23), the proof is complete. □ 

Now we turn to the multidimensional CLT for random vectors, 

S s {X 1 W n ) = (MAu 1 (X 1 W n )),...,iy d (A Ur (X,W n ))) T 1 n e N, (28) 

where u= (iti, . . . , u,-) T £ M. r . 

Theorem 2. Let X = {X(t),tE M. d } be a random field satisfying condition (A). Then, 
for each u— [u\, . . . , u r ) T £ R r and any VH-growing sequence (W n ) T igN of subsets ofM. d , 
one has 

VdiW^^iS^X^-vaiW^Piu^AVz-N^iu)) asn^oc, (29) 

where 

P(u) = (P(X(0) > ui), . . . , P(X(0) > u r )) T 
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and £(u) = (c; TO (u))j" =1 is an (r x r) -matrix having the elements 

ai m (u)= [ cov(I{X(0)>ui},I{X(t)>u m })dt. (30) 



Proof. Observe that the convergence of all r 2 integrals in (30) is proved in the same 
way as that of the integral representing cr 2 (u) in the one-dimensional case. The result 
follows by using the Cramer-Wold device. We omit further details that are quite similar 
to those in the proof of Theorem 1. □ 

The last theorem entails: 

Corollary 1. Let X = {X(t),t G R d } be a random field satisfying condition (A). Assume 
that X(w) is non- degenerate for some u G W . Then, for this u and any sequence (W n )ne?$ 
of VH-growing subsets ofM. d , one has 

^(V^r 1/2 £(«r 1/2 (^(X, W n ) - v d {W n )P{u)) AV~ 7V(0, 1) as n^oo; 
here I denotes the unit (r x r) -matrix. 

3.2. Shot noise random fields 

We verify the conditions of Theorem 1 for shot noise random fields. These fields appear 
naturally in the theory of disordered structures. Let 2?(R d ) (resp., Bo(M. d )) be the family 
of all (bounded) Borel sets in M. d . A shot noise random field X = {X(t),t £ M. d } is defined 
by the relation 

*(*) = £&¥>(*- a*)> 

where {£»} is a family of i.i.d. non-negative random variables and {x{\ is a homogeneous 
Poisson point process in M. d with intensity A <E (0,oo), that is, {x{\ is the support set 
of a random Poisson counting measure {Nb, B g B(M. d )}, where Nb = #{i: x, £ B} has 
the following properties: 

(i) Nb i: Nb 2 :--- are independent for pairwisc disjoint Bi,B2, . ■ . G Bo(M. d ), 
(ii) N B - Pois(A^ d (B)) for all B G B (R d ). 

Suppose that {£j}, N/.\ arc independent, E^f < od and ( y 9:M c '^M + is a Borel function. 
For the shot-noise field X introduced above, we impose the condition: 
(C) X(0) has a bounded density and for a function ip bounded and uniformly contin- 
uous on IR d , 

¥>(*)< W>(ll*ll2) = 0(||t||j a ) as Ht|| 2 ^oo, (31) 

where a > 3d and ipo ■ R+ — > M+ ■ 
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Proposition 1. The statement of Theorem 1 holds for a random field X satisfying 
condition (G). 

Proof. By [7], Theorem 1.3.8, X is associated and hence quasi-associated. Moreover, it 
is strictly stationary with covariance function given, for example, in [7], Theorem 2.3.6. 
The continuity of the covariance function follows from the inequality 

| cov(X(0),X(s)) - cov(X(0),X(t))| < AE£ 2 sup \<p(t - y) - <p(s - y)\ f \<p(y)\dy 

y<ER d JR d 

and the uniform continuity of ip. Corollary 2.3.7 of [7] yields the desired bound for the 
covariance function in condition (A) . The proof is complete. □ 

Note that the characteristic function of X(0), provided by [7], Lemma 1.3.7, is inte- 
grablc if 

explx (<pt:(sip(t)) - 1) dt\ ds<oo. (32) 

Thus, (32) guarantees the existence of the bounded density of X(Q). 

Condition (32) can be easily verified in a number of special cases; for instance, if 
£i = const > a.s. and ip(t) = aexp{ — 6||i|| 2} or ip(t) = amin{l, ||i|| J } with a, b > 0. 

3.3. Gaussian random fields 

In contrast to Lemma 1, we obtain a sharper estimate for the covariance of indicator func- 
tions in the Gaussian case. Our result extends formula (2.7.1) of [12]. Let $ and V? stand 
for the cumulative distribution function and the tail distribution function of a standard 
Gaussian random variable, respectively. 

Lemma 2. Let (U, V) T be a Gaussian random vector in R 2 such that U ~ J\f(a,T 2 ), 
V ~ 7V(a,T 2 ), where a£l, r > and correlation coefficient corr(U,V) = p. Then, for 
any B^et and p € (— 1, 1), the following eguality holds: 



co\i{I{U>u},l{V>v}) 

1 f 1 f (u-a) 2 -2r(u-a)(v-a) + (v-a)' 2 . 

ex P KZTT, 15* t dr - 



(33) 

2nJ vT^P^l 2r2(l-r 2 ) 

In particular, for u~v, one has 

„v(. {I />„ } ,.{V>»})^f 7I ^ s «p{-^}d, 

Moreover, for any u,v € R and p G [— 1, 1], one has the ineguality 

|cov(I{U">it},lI{F>v})|<|p|/4. (34) 
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Proof. Using the transformation ih-(j;- o,)/t, x £ R, we can assume w.l.g. that U ~ 
7V(0,1) and V~jV(0,l). Let p£ (-1,1). The probability density 

1 f x 2 -2pa;y + y 2 

fu,v(x,y) = , = exp' 



2n^7 I 2 (!-P 2 ) 

of the bivariate Gaussian random variable (£/, V) is invariant under the transformation 
£i-> —x and yH> — y, (x 7 y) T £ R 2 . Therefore, 

cov(I{L7 > u},I{V> u}) = cov(I{£7 < -u},I{V < -u}), u,v£R. 

It is well known (see, e.g., [10], formulae (21.12.5) and (21.12.6)) that 

~ $( fe +i)(y ) $(*+i)(V ) , 

/^(z,y) = £- ^f ^V, ^£R, 

fe=0 

where $( fc ) (x) = d k <&(x)/dx k and, for any u, v £ R, 
/ /c/.v(a;,2/)d.Tdy = ^ 

-OO «/ — OO (. n 

Hence, for each u, u € R, 



00 - 7 - 00 fc=0 



$(*)( u )$(*)(«) k 



/ f uy {x 1 y)Axdy = Y J 

OO J — OO 7. n 



$(k) (_«)$(*) (_ t 



oo-z-oo fe=Q 

fe=l 

*/ x*/ ^ rA$( fe + 1 )(- U )cf>( fc + 1 )(- U ) . 
= $(-«)$(— v) + / 2^ jp -r k dr 

*' fc=o 

= $(-u)$(-v) + / /[/(r),V(r)(w,v)dr, 

Jo 

where (C/(r),T^(r)) T is a centered bivariate Gaussian vector with EU(r) 2 = EV(r) 2 = 1 
and cov([/(r), V(r)) = r. Consequently, we get 

cov(I{[/< -u},l{V< -v}) = / fu(r),v(r)(u,v)dr 





1 /•" 1 



zC -(u 2 -2ruv+v 2 )/(2(l-r 2 )) ^ 



2tt7 VT^T 2 

Passing to random variables U and V with arbitrary mean a and variance r 2 > gives 
the formula (33). 
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To prove inequality (34) for p e (—1,1), write 

\cov{l{U<-u},I{V<-v})\<^J _L_dr<i-arcsin|p| 

and notice that arcsin|/9| < n\p\/2. 

The case \p\ = 1 is trivial, as \P(AC\B) - P(A)P(B)\ < 1/4 for any A, Be J 7 . □ 

The following result generalizes the corresponding one established in [12] (see page 80), 
where the isotropy of the Gaussian random field was assumed. A central limit theorem 
for nonlinear transformations of a homogeneous Gaussian random field was used there. 

Theorem 3. Let X = {X(t),t e K d } be a Gaussian stationary random field satisfying 
condition (B) and X(0) ~ J\f(a,T 2 ). Then, for each u£l and any sequence of VH- 
growing sets W„ C M. d , one has 

v d (A u (X, W n )) - v d {W n )*{{u - a) /t) ± , 

i — , TT . N -^ i u ^ j\ y\J , a \u)) 

\JVd{W n ) 

as n — > oo. The variance cr 2 (u) introduced in (11) can be written in the following form: 
* 2 M = ±-[ f (t) -^exp(-^-^-W (35) 



ZnJwJo VT^ "\ r 2 (l+r) 
where p(t) = corr(X (0) , X (t)) . In particular. 



a 2 (a) = — / arcsin(p(i)) dt. 

27t JRi 



Proof. For the Gaussian field X, we have P(X(0) > u) = ty((u — a)/r). Now we ap- 
ply the upper bound (34) to obtain \cov{X(0),X(t))\ instead of | cov(X(0),X(i))| 1 / 3 in 
the estimates used in the proof of Theorem 1. This leads to the hypothesis that a > d 
in (9) , whereas in condition (A) we assumed a > 3d. Note that Gaussian fields are quasi- 
associated [20]. 

Finally, we express <t 2 (m) (see (11)) in terms of the covariance function of X as in the 
proof of Lemma 2, and this yields (35). The proof is complete. □ 

Theorem 4. Let X = {X(t),t e M. d } be a random field satisfying condition (B) and 
X(0) ~ Af(a, t 2 ). Then, for each u = (ui, . . . ,u r ) T el'' and any sequence (W n )neN of 
VH-growing subsets ofW 1 , one has 

v d {W n )- x ' 2 {S u {X,W n )-v d {W n )*(u))AVz~N{Q,Z{ti)) asn^oo. (36) 

Here, \&(u) = {^{{u\ — a)/r), . . . , $((u r — a)/r)) T and £(w) = (c/ m (M))[ m=1 is a matrix 
having the elements 

I r rp(t) 
o-i m (u) = — / g(r)drdt, (37) 

2n Jw Jo 
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where 

1 f (ui -a) 2 -2r(ui -a)(u m -a) + (u m - af 

9{r) = 7T = ^ CXP { 2^1=^ 

and p(t) = corr(X(0),X(t)). //£({?) is non- degenerate, we obtain by virtue of (36) 

^(W„)- 1 / 2 S(u)- 1 / 2 (^(X, W n ) - v d {W n )^{u)) A AA(0, 1), n -> oo, 

where I is the unit (r x r) -matrix. 

Proof. Employing Lemma 2, one can repeat the reasoning proving Theorem 2. Clearly, 
P(X(0) > u{) = *&((ui — o)/t). I = 1, . . .,r. The matrix elements ai m (u) for l,m= 1, . . .,r 
can be calculated by way of (33). □ 

Formulae in the isotropic case 

For the isotropic case, we use the change of variables (passing from t= (ti,... ,^d) T to 
spherical coordinates) in integrals (35) and (37) to obtain the following statement. 

Corollary 2. Let the random field X — {X(t),t € R d } satisfying the conditions of The- 
orem 3 be isotropic (d>2). Then 

a , . duj d f°° d , T^ 1 f (u-a) 2 1 , , 

o- 2 t u ) = —± v d - 1 exp^ — Vt- — ^-Urdu, 

1 ^ 27t 7 i VI372 *\ r 2 (l + r)J 

where p{v) = corr(X(0),X(t)) if\t\—v. Foru — a, one has 

o- 2 (a) = — — / w d_1 arcsin(p(-u))dt;. 
2n Jo 

In the multivariate case, (37) can be written as follows: 

^m(«) = -^7 / w / g(r)drdw 
/or Z, 77i = 1, . . . , r. 



27T Jo Jo 



4. Statistical version of the CLT 

Now we provide a statistical version of the CLT involving random self-normalization. Let 
r £ N be the number of levels to observe. 

Theorem 5. Let X = {X(t),t G R d } be a random field satisfying condition (A). Let Uk € 
M, fc = l,...,r and (W n )„ e n be a sequence of VH-growing sets. Furthermore, let C n — 
(c ra Zm)[ m -i be statistical estimates for non-degenerate asymptotic covariance matrix E 
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with elements <x\ m given by (30). Assume that c n i m — > o~i m as n — > oo for any l,m = 
l,...,r, where — > denotes convergence in probability. Then 

C- l ' 2 v d {W n y l l 2 {S{W n ) - v d {W n )P{u)) A N{0,1) asn^oo. 

Proof. It suffices to use Theorem 2 and elementary properties of the convergence in 
probability and in law for random vectors. □ 

One feasible estimator for the asymptotic covariance matrix £ that arose in the multi- 
variate CLT, see Theorem 2, can be called a subwindow estimator [18] and is constructed 
as follows. Let (V n ) n ^ and (W^)n€N De sequences of VH-growing sets (not necessar- 
ily rectangles) such that V n C W n , n> 1. Consider N(n) subwindows V n ,i,..., V n ^( n ), 
where (N(n)) n£ jq is an increasing sequence of integers with lim ra _> oo N(n) = oo, and 
Vn,j =V n ® {ftuj} are subwindows that are translated by certain vectors h n j £ M. d , 
j = 1, . . . , N(n). Assume that U 7 =i ^»,j — W« f° r each neN and there exists some 
r > such that 

V nJ n V nii C 9Kj © B r (0) for i, j G {1, ■ • ■ , N(n)} with i ^ j. 

Denote by 

£$ = -^tt / HX(t) > u k } dt, j = l,..., N(n), 
Vd{V n ) J Vn . 

the estimator of /j,k = P(X(0) > Uk) based on observations within V n j, and by 

JV(n) 

Jink = Wn) S#£k' «GN, * = l,...,r, 

the average of these estimators. After all, we define the estimator £„ = (oViim)[ m=1 for 
the covariance matrix £. Set 

i/ rv ) N{n) 

&nlm = JV(n) - 1 ^ ^ ~ A™^^"™ ~ A"™)' ^ 38 ^ 

We recall the following result. 

Theorem 6 ([18], Theorem 3). Let X = {A(i),iG IR d } be a strictly stationary random 
field such that 

I (2 2) 

/ Rm ( a; >2/^)l dxdydz<oo, f,m = l,...,r, (39) 

where the fourth-order cumulant function 

(2 2\ 

c \ m (x,V,Z) = E([^(0) - in][Z m (x) - ^ m ][Z ; (y) - m][Z m (z) - /x m ]) 

- cov Zm (ir) cov im (z - y) - cov u (y) cov mm (x - z) - cov im (z) cov mi (a; - y) 
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and cov; m (i) = cov(I{X(0) > ui},I{X(t) > u m })> l,rn= l,...,r. Then £ introduced 
in (38) is mean-square consistent. 



Relation (39) holds for a random field X with finite dependence range. In this case, 
the estimator £ n is mean-square consistent. Among other estimators for the asymptotic 
covariance matrix, there are two worth mentioning. One estimator that, under certain 
assumptions, meets the conditions of Theorem 5 is introduced in [6, 8] and involves local 
averaging. A major disadvantage is tedious calculation in the case of a large observation 
window. The same problem arises for an estimator based on the covariance function 
estimation for the underlying random field; see [18] (cf. [12], Chapter 4). 



5. Discussion 

A very important issue for applications of the estimator E n is the choice of an appropriate 
size of the (e.g., rectangular) subwindow V n . The subwindow size is related to both 
the covariance structure of the considered random field and the size of the observation 
window. We will discuss these problems while considering a simple example. The data 
used consist of 100 mutually independent realizations of stationary and centered Gaussian 
random field X with covariance function 

cov(*(0),*(t)) = (l- ^ + ^) 1 {ll*lla e [0,0]}, t e K 2 , 

for some a > according to the spherical covariance model (sec [21], page 244), which 
is often applied in geostatistics. The correlation range a in our simulation study has to 
be small enough in comparison with the size of the observation window to make the 
CLT argument work. Here, we take, for example, a = 10. The fields are simulated in the 
observation window W = [0, 2000) x [0, 2000) on the grid with mesh size one. That means 
every realization provides 4 million data points. To generate level sets, we consider the 
thresholds u\ = —1.0, «2 = 0.0 and U3 = 1.0. Then 

'4.6432 

I 5.9938 10.5564 
2.7962 5.9938 4.6432, 

An appropriate subwindow size can be found focusing only on the threshold «2 = 0.0, 
since for other threshold values the obtained results differ from this one only slightly. The 
estimator provides the best result for £ as the edge length of the rectangular subwindow 
equals 15. In general, the optimal choice of this length is an open non-trivial problem. Af- 
ter this preliminary step, we are able to apply the subwindow estimator to the simulated 
data. The following two matrices show averaged estimation results for £ by means of £. 
On the left-hand side, the averaged value of each estimated matrix element is computed 
out of 100 samples. On the right-hand side, the mean error to the theoretical value is 
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provided. 

x 100 / 4.6556 \ / 0.27% 

-— ^ S fc = 5.9710 10.5524 , ME = -0.38% -0.04% 

100 fc=i \ 2.8156 5.9934 4.6762/ \ 0.69% w 0.00% 0.71%, 

It would be interesting to propose a statistical hypothesis test based on the established 
statistical version of the CLT in order to apply it to data concerning the paper production. 
We will deal with this topic in a separate paper. 

6. Open problems 

The research area of limit theorems for level sets of random surfaces still offers an abun- 
dance of open problems. Let us mention just a few. It would be desirable to prove limit 
theorems for joint distributions of various surface characteristics of different classes of 
random fields. For instance, one could consider stable fields. Further on, one can study 
random fields possessing more strong dependence structure: for example, satisfying con- 
dition (A) with a < d. In this case, the normalizing factors have to be changed and the 
limiting distributions can be non-Gaussian. Certain results for problems of this type can 
be found in [12, 15]. One could also prove a functional limit theorem for an innumerable 
set of thresholds. As our main result could also be called the CLT for the first Minkowski 
functional, it might be of interest to prove limit theorems involving other Minkowski 
functionals for level sets such as the boundary length or the Euler characteristics. It is 
worth mentioning that, for a stationary two-dimensional Gaussian field, this has already 
been done for the second Minkowski functional in [14]. 
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